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Abstract 



£5 ' We carefully investigate the two fundamental assumptions in the Stillinger- 

O ' Weber analysis of the inherent structures (IS's) in the energy landscape 

and come to conclude that they cannot be validated. This explains some 
of the conflicting results between their conclusions and some recent rig- 
orous and exact results. Our analysis shows that basin free energies, and 
^ ■ not IS's, are useful for understanding glasses. 

c«v 

' It is well known that most supercooled liquids (SCL) become viscous when 

\Q , their configurational entropy S(T) becomes negligible as they are cooled, 

provided the corresponding crystal (CR) is not allowed to nucleate. The flow 
practically ceases over a period constrained by experimental limits, the viscosity 
becomes very large, and the viscous fluid eventually becomes an amorphous solid 
or glass. Our current understanding of glassy behavior is still far from complete, 
even after many decades of continuous investigation. In order to better under- 
stand the flow properties of viscous fluids, Goldstein proposed the potential 
I ■ energy landscape picture using classical statistical mechanical canonical ensem- 

ble 2 , in which the energy barriers control the flow at low temperatures. The 
discussion was mostly qualitative, but provided an interesting and sufficiently 
£j \ tractable scheme and included some quantitative predictions in the SCL and 

the glassy states. Stillinger and Weber (SW) later revived this picture, carried 
out an analysis in terms of basins, and concluded that their minima, called the 
■ inherent structures (IS), play a pivotal role in the thermodynamics of viscous 

H | fluids at low temperatures 0] . The IS-entropy Sis (T) [not to be confused 

with S(T)] of the IS's vanishes at some T = T sw (called T K by SW, but we will 
reserve Ik where S(T) = pQ) so that the system gets trapped into a single 
basin at Tsw so that S(T) becomes the basin configurational entropy Sh{T) > 
below Tsw > ^kOQ- The SW analysis has given rise to a considerable amount 
of literature in recent years; for a partial list, see [S]. 
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One of the most tantalizing consequences of the SW analysis is the theoret- 
ical conclusion Tsw = drawn by Stillinger 4 . Consequently, it is not possible 
to have a configurational entropy crisis (S(T) < 0) below a positive tempera- 
ture Xk PP for SCL. For CR, the claim implies that just above Tsw = 0, CR 
probes many basins (Sis > 0), not all of which have to be close in the configu- 
ration space, and its heat capacity is not due to pure vibrations within a single 
basin. This is hard to understand in view of the tremendous success of the 
Debye model. Stillinger's argument also does not permit any SCL spinodal at 
a positive temperature. However, the latter has been observed in exact calcula- 
tions for finite-length polymers jHj and in a binary mixture [7]. The theoretical 
conclusion Tsw = 0] is contradicted by numerous numerical evidence of 
a positive Tsw within the IS picture. A recent rigorous statistical mechanical 
proof by Gujrati 3] also contradicts this claim. He has shown that under the 
assumption that the ideal glass has a higher energy (E = Ek) than the corre- 
sponding crystal (E = Eq < Ek) at absolute zero, S(T) of the stationary SCL 
(obtained under infinitely slow cooling of the disordered equilibrium liquid EL) 
must necessarily vanish at Tk > 0. At Tk, SCL has its configurational energy 
Ek- Two independent proofs are given jS], and the conclusions are also substan- 
tiated by two exact model calculations, one of which is not mean-field. Other 
exact calculations for the abstract random energy model [5], and for long poly- 
mers jHlEl a l so support the conclusions by Gujrati. In addition, a recent exact 
solution by Semerianov and GuiratifTT] of a dimer model, a prototype model of 
molecular liquids, also exhibits configurational entropy crisis in its stationary 
metastable state below a positive temperature. Corsi ^2] has also observed 
positive Kauzmann temperatures in exact calculations for small particles oc- 
cupying four and fives lattice sites each. The important point to note is that 
none of these exact calculations utilizes the energy landscape. However, if the 
latter has any validity, its consequences must be in accordance with the exact 
calculations and the rigorous analysis, which is certainly not the case. Thus, 
we need to reexamine the SW analysis to clarify the conflict. We also provide 
an alternative analysis of the landscape, which is consistent with the rigorous 
analysis. 

The canonical ensemble free energy -Fdis(T) of the disordered EL is continued 
analytically [Sj below the melting temperature Tm to give the SCL free energy. 
As a mathematical continuation, the resulting SCL does not have to satisfy 
the reality condition S(T) > 0. The continuation, in principle, can stop in a 
spinodal singularity at a positive temperature [HIE]. Here, we are only interested 
in the case when the SCL free energy can be mathematically continued all the 
way down to T = without encountering any singularity. It is easy to show 
|S] that, at T = 0, the SCL and CR free energies are identical (F = E ), 
provided TS(T) — > 0. They are again equal at Tm because of which SCL is 
forced to exhibit the entropy crisis (S(T) < 0) below T K > ED EH, 
where it is discarded and replaced by an ideal glass phase to satisfy the reality 
condition. The ideal glass has E — Ek for T < Tk- At T = 0, the configuration 
corresponding to Ek must represent a potential energy minimum, provided we 
neglect surface effects. At a positive temperature Ts > Tk, SCL gets confined 
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into the basin whose IS is at . Thermal fluctuations will still allow SCL to visit 
other basins, which are rare if T is small. We will not investigate fluctuations 
and restrict ourselves only to the average behavior here. 

The following two observations jS] are going to be relevant below. 

Gl. The zero of the temperature scale is set by the global potential energy 
minimum Eq. 

G2. The ideal glass configuration at -Ek is a local minimum of the potential 
energy (neglecting surface effects), which SCL approaches at Tk > 0. Since SCL 
does not physically exist below Tk (due to negative entropy), the ideal glass is 
put in by hand for T < Tk to avoid the entropy crisis. It does not emerge 
directly in the statistical mechanical description. 

In continuum classical statistical mechanics, the dimensionless total partition 
function (PF) Zt of a system of N identical particles i — 1, 2, N in a given 
volume V, with the Hamiltonian £/r({i"i, Pi}) = K({pi})+E({ri}), K,E being 
the kinetic and the configurational (i.e., the potential) energy, is a product of 
two independent integrals 



Here d N {p},d N {r} represent integrations with respect to momenta and posi- 
tions Pi,Ti of the particles, and (3 the inverse temperature 1/T in the units of 
the Boltzmann constant fee- The second integral (the prime implying integra- 
tion over distinct configurations of the particles) is called the configurational 
PF, to be denoted by Z. 

The momentum integration in Z^ can be expressed in terms of WKE(P)dP 
= C3jyP 3N ~ 1 dP/h 3N related to the 37V-dimensional momentum space volume 
within the spherical shells of radii P, and P + dP [C d = dir d ^/T(d/2 + 1)]. The 
translational entropy due to the translational degrees of freedom is given by 
Ske(T) = lnW KE (P), K = P 2 /2m = 3NT/2. In the thermodynamic limit 
N -> oo, we find that S K e(T) = (3iV/2)[l + InT + ln(27rm//i 2 )], and has 
the same value at a given temperature for all classical systems, regardless of 
their configurational energy. Thus, in general, the entropy due to the config- 
urational degrees of freedom can be always obtained by subtracting Ske(T) 
from St(T)H3 E] : S(T) = S T (T) - S KE (T), where S T (T) is the total en- 
tropy in the canonical ensemble; see (Q. For an ideal gas, Z = V N /N\, so 
that S = N\n(Ve/N), which no longer depends on T. For T < eh 2 /2Trm, 
S K e(T) < 0, or for V/N < 1/e, S(T) < 0. 

The problem of negative entropy is well known in classical statistical me- 
chanics, and requires quantum statistical mechanics for its resolution. Unfortu- 
nately, it is not possible to solve a quantum statistical mechanical model exactly 
at present. Thus, care must be exercised when drawing conclusions based on the 
sign of the entropy. Since Tk is signaled by a negative entropy, it is crucial to 
have a formalism in which the entropy is never negative for realizable states 
in Nature EI EH]- The simplest way to achieve this is to discretize either 
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the phase space by using cells of size h 3N or the real space by using a lattice. 
Thus, in the following, we take it for granted that such a discretization has been 
carried out. We closely follow ^ who factors out the kinetic energy part and 
only uses Z to introduce the potential landscape picture [21 EH, and consider only 
the configurational entropy S(T) pQ, the canonical PF Z(T), and the configu- 
rational free energy is F(T) = —Th\Z(T) in the following. We rewrite Z(T) 
as follows. Let W(E) > 1 (so that S(E) = ]nW(E) > 00 E3| for physically 
realizable states) denote the number of the configurations of potential energy 
E. Then Z(T) can be rewritten as a sum over E : Z(T) = ^ B e^ B lf(£). 
Replacing W(E) in Z(T) by the number of disordered configurations W<a s {E), 
we obtain the PF Z dis (T) so that ^ d is(T) = -TlnZ dis (T) is the free energy of 
the disordered phase EL and its extension SCL [3J El Ej • In the following, we 
will use Z(T) to represent both PF's, which should cause no confusion since the 
context will be clear. 

SW Picture. The potential energy landscape is a union of disjoint basins. 
A basin is indexed by j, and characterized by its minimum and maximum en- 
ergies Ej, and Ej >max so that it does not exist outside this energy range AjE = 
(Ej,Ej <max ). Let Wj(E) (E e AjE) represent the number of distinct states of 
energy E in the j-th basin. We introduce the shifted PF 

Zj (T) ee YWiiE)e- 0{E - Ei) Sea,e (2) 

E 

of the j-th basin. Here, Se,AjE = 1 if E € AjE, and if E £ AjE. We group 
basins, indexed by j(X), into inherent structure classes (ISC) X x , indexed by 
A, so that all IS's in a class have the same energy E = E\. The basin in a 
class do not have to be close in the configuration space. Let Ni$(E\) be the 
number of basins in 2a, and S\s(E x ) = lnNis(E\). Let Z\ = J2jej{\) Z A^) an< ^ 
z\ ee Z\/Nis{E\) denote the shifted and the average shifted I\-P¥, so that 

Z(T) EE ]T e-^Z x EE e-^ +Sls(E ^z x . (3) 

A A 

Stillinger and Weber [3JE] and various authors in [S] assume that z\ is an explicit 
function of E\, and T of the form [SW denote quantities that are specific to the 
SW-approach] 

zf W \E^T)^e-^ SV,) ^ T \ (4) 

They also replace the first sum over the discrete index A in Q by a sum over the 
almost continuous variable E\ so that a general summand can be characterized 
by E x using F( SW \E X ,T) ee E x + /< sw ) (E X ,T) - TS 1S {E X ), whose minimum 
with respect to E x at fixed T determines Z(T) for a macroscopic system. This 
minimum term corresponds to that particular value E x = E x at which 

(8S IS (E x )/dE x )z x = (3[l + (df( sw \E Xl T)/dE x )_ }. (5) 

\ / E\ 

The equilibrium free energy, and the IS-entropy are given by F^ SW \t) ee 
F^iE^T), and S^ W) {T) = Sis(E x ), respectively. We expect 3^ (T) 
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to vanish at some low T = Tsw and increase monotonically with T at least at 
low T. For T < Tswj the system is trapped in a single basin of energy minimum 
E x = E K for SCL (E for CR) so that at or below T sw , E x sticks at E K for SCL 
(E for CR). Stillinger ^ has argued that satisfying © below Tsw is inconsis- 
tent with E\ sticking at either Ek (or E ). Thus, he concludes that Tsw = 0, 
a consequence of which is the observation that there cannot be a positive Tk- 
This conclusion is based on the following two fundamental assumptions in the 
SW approach: 

SW1: The shifted ISC free energy /( sw )(£: a ,T)-TS*is(Ta) in T< sw ) (E x , T) 
is an explicit function of E\ , so that the summation over A can be replaced by 
that over E\. 

SW2: The minimization condition (j^J also holds below Tsw, where 

(). 

We now demonstrate that neither assumption can be substantiated. 

Current Analysis. We observe that W(E) is a sum over various basins: 
W{E) = Y^j Wj(E)5E,AjE- We group all basins having the same minimum at 
E = E\ together into X x . It should be noted that Zj (T) of all basins in T x need 
not be the identical in value. 

1. Evaluating zj(T). We now prove that Zj cannot depend explicitly 
on the basin energy minimum Ej , see (|2J , though it most certainly depends on 
the shape of the basin, i.e. on j. For example, the curvature of the basin at 
its minimum and not its value of Ej determines the vibrational frequencies and 
the free energy fj(T) = — Tlnzj in the harmonic approximation. The latter is 
measured with respect to Ej, so is independent of Ej. To be sure, let us shift 
all energies E — > E = E — C by some constant C in J2J). The number Wj(E) 
of states, all having the same energy E, remains unchanged under the shift by 
C. Thus, Wj(E) -> W'j{E') = W 3 {E). Thus, z 3 transforms under the shift as 

Zj(T) — > J2e' Wj(E )e~ f3l - E ~ E ^8 E ' A . E > for any arbitrary C. Comparing with 
J2J), we conclude that Zj has not changed. Consequently, it does not depend on 
the shift C, including C — Ej. 

Since Wj(E) > 1, Zj(T) is a sum of positive terms. Hence, for a macro- 
scopic system, Zj is determined by the maximum summand in @ correspond- 
ing to E — Ej G AjT', and the corresponding heat capacity is non- negative. 
(Both observations remain valid even if Wj(E) > 0, a common occurrence in 
SCL continuation |H1 151 ITU1 HT1 112] .1 For E = ~Ej. we have 

(dS 3 (E)/dE) Ej =(3, % G AjE or T G A 3 -T, (6) 

where AjT is the temperature range (Tj-, Tj, max ), so that the equilibrium energy 
Ej (T) for the basin lies in the range AjE. According to G2, Tj- is strictly positive 
for Ej = Ek > Eq. Thus, we assert that Tj is not necessarily zero The 
basin free energy is 

fj (T) = -T In Zj = Ej - Ej - TS, (Ej), T G AjT. (7) 
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The energy landscape is topologically very complex, with various basins very 
different from each other, even if they have their minima at the same energy. 
Thus, Ej(T) in different basins at the same temperature T (provided T £ AjT 
for these basins) will be usually different. There is no requirement that they 
be the same. Moreover, even if the free energies fj(T) of two or more basins 
happen to be the same at some temperature T, they need not remain equal at 
other temperatures. 

2. Evaluating Z X {T). The proper form or the ISC PF's Z X {T) or z x in 
T\, see ©, is 

Z x ^z x e s ^^ Z ^A ' • (8) 

Due to the delta term, similar in definition to Se,AjE, the sum in JSJ contains 
only those basins in T x that exist at T in the sense that its range contains 
T : AjT 3 T. From now onward, we only consider those basins that exist 
in this sense; hence, we will not explicitly exhibit the delta term anymore. 
We now classify each existing basin in (JSJ according to its free energy /. Let 
N x (f) denote the number of basins of free energy / that exist at a given T in 
I x . Since / is a function of T, S x (f) = lnAT^(/) also changes with T. For a 
macroscopic system at a given fixed T, Z x is dominated by the basins in I x 
for which F x (f, T) = f — TS X is minimum as a function of / at fixed T. The 
resulting entropy and the free energy at the minimum (/ = f x ) are denoted by 
S X (T) = S x (f x ), and F X (T) = f x — TS X , respectively, where f x is determined 

by 

(dS x /df) Jx =(3, T>T^ X) . (9) 

which looks similar to © but very different from J3J). Here, lg ^ is the tem- 
perature at which S x = 0, so that only one basin exists below it (and above 
its lowest temperature Tj) in the above sense. Because of this, the issue of 
F\-niinimization for T < does not arise as there is only one member in X x 
(S x = 0), so that F x — > f x ; thus, the minimization of Fa is already ensured by 
(0. We should contrast @ with JSJ. Stillinger puts no restriction on the ap- 
plicable temperature range in the latter; see SW2. However, it is clear from our 
discussion that (jSJ) cannot apply below Tgw, which invalidates the second fun- 
damental assumption SW2. Consequently, Stillinger's argument that Tsw = 
has no validity 0]. 

3. Evaluating Z(T). Using the evaluated Z x in Z, we find 

Z = Y,e-ft E *+ F *™. (10) 

A 

We now make the following important observation. As shown above, f x is 
independent of E x ; thus, S X (T) and F X (T) are independent of E x , although 
they most certainly depend on the ISC X x . There are various aspects of the 
basins such as the IS-curvature that determine F x , but E x is not one of them. 
This has the following very important consequence. From JHJ, we note the 
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identity f x = F x + TSi S (E x ), where f x (E x ,T) = —T In z x [ compare with 
f@W)(E x ,T) in Q]. The explicit dependence of f x {E x ,T) [or /( sw ) (E x , T)] 
on i?> must be trivial and due to S\s(E x ) because of the independence of F x 
[or f( sw \E x , T) — TS\${E X )] on E x . Accordingly, the summation over A cannot 
be replaced by a summation over E x . This disproves SW1. 

We now deal with the summation over A in i|l(J|l in a standard manner. Let 
denote the number of ISC's at a given T with the same F =E X + F X (T). 
The PF Z is dominated by the ISC's of free energy T for which T — TS is 
minimum over J-; here S(J-) = In M{T). The condition for this minimum at T 
is 

(dS/dT) T = /3, T >T, (11) 

as expected. Here, T is the temperature at which S = S(J-) = 0. Hence, we 
finally conclude that the final free energy is given by 

F = -TlnZ =7 -TS, T > T b , (12) 

where the significance of Tb will become clear in a moment. It should be obvious 
at this point that the dominant contribution in l|10|) mixes ISC's with different 
E x and Tg . For T < T, the system is confined to a single ISC corresponding 
a particular value A = a. The minimization of T — TS is no longer an issue for 
T < T , as there is only one ISC a to consider and the minimization of T = F a 
is already ensured by @. 

It is easy to see that the configurational entropy is 

S(T) = S + S a + S h . (13) 

As the temperature is lowered, regardless of whether we consider CR or SCL, we 
must first encounter the case S = at T, so that the system is confined to the 
single ISC a as discussed above. At a lower temperature Ts = Tg^ <T,S a = 0, 
and the system is trapped in a single basin j =b. (Compare Ts with Tg\y.) The 
basin j =b has its lowest allowed temperature T = Tb < Ts, where the basin 
entropy S*b = 0. The configurational entropy vanishes when all three components 
in (^21 vanish. This obviously happens at T = T b < T s . For CR, T b = 0, and 
for SCL, Tb = Tk > 0, as proven recently|S]. We now restrict our discussion 
to SCL and consider Tb < T < Ts, so that S c — and the system is trapped 
in the basin j =b whose IS is at Ek and Tb = Tk. In this basin, © remains 
satisfied until T > Tk, but ceases to work below Tk, if we insist that E < Ek is 
not allowed. However, if we do allow E < Ek, which requires a mathematical 
continuation of SCL free energy function below Tk (which is not equivalent to 
the continuation of the SCL state itself as it does not exist below Ek), we can 
continue to impose ©. This will produce a negative entropy continuation of 
the SCL-entropy down to absolute zero, which is what all analytical calculations 
show in which the SCL free energy is obtained by continuing the disordered 
phase free energy below Tm- However, demanding that 10 continue to operate 
below Tk, and at the same time demanding that E remains fixed at E=Ek is 
not legal. Similar arguments apply to the consideration of S and Ts. We cannot 
demand to remain operative below Ts; it must be replaced by ijfijl. 
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It should be noted that the general landscape picture itself is not capable 
of providing any information about whether Tk > or not. For this, we must 
turn to other approaches like the one developed by Gujrati|8] or to some model 
landscapes. Indeed, model landscapes can be easily constructed in which the 
temperature associated with the minimum of a basin need not vanish |13) . 

In summary, we have shown that when the landscape picture is carefully 
developed, there is no conflict with the exact and rigorous results about Tk. In 
the process, we have also corrected some of the flaws in the SW analysis. As the 
current analysis deals with the free energies and not the energy minima, IS's 
play no useful role except possibly at very low temperatures below Tg. 

It is our pleasure to thank Andrea Corsi for his comments on the work. 
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